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Non-ionized media subject to strong fields can become locally ionized by penetration of finger- 
shaped streamers. We study negative streamers between planar electrodes in a simple deterministic 
continuum approximation. We observe that for sufficiently large fields, the streamer tip can split. 
This happens close to the limit of "ideal conductivity". Qualitatively the tip splitting is due to a 
Laplacian instability quite like in viscous fingering. For future quantitative analytical progress, our 
stability analysis of planar fronts identifies the screening length as a regularization mechanism. 



Streamers commonly appear in dielectric breakdown 
when a sufficiently high voltage is suddenly applied to a 
medium with low or vanishing conductivity. They consist 
of extending fingers of ionized matter and are ubiquitous 
in nature and technology ||l|,|l. The degree of ionization 
inside a streamer is low, hence thermal or convection ef- 
fects are negligible. However, streamers are nonlinear 
phenomena due to the space charges inside the ionized 
body that modify the externally applied electric field. 
While in many applications, streamers by a strongly non- 
uniform background electric field are forced to propa- 
gate towards the cathode through complex mixtures of 
gases (Ij-Q, we here investigate the basic phenomenon 
of the primary anode-directed streamer in a simple non- 
attaching and non-ionized gas and in a uniform back- 
ground field as in the pioneering experiments of Raether 
In previous theoretical work, it is implicitly assumed 
that streamers in a uniform background field propagate 



in a stationary manner . This view seems to be sup- 
ported by previous simulations [p|jic|]. 

In this paper we present the first numerical evidence 
that anode directed (or negative) streamers do branch 
even in a uniform background field and without initial 
background ionization in the minimal fully determin- 
istic "fluid model" if the field is sufficiently 
strong. We argue that this happens when the streamer 
approaches what we suggest to call the Lozansky-Firsov 
limit of "ideal conductivity" The streamer then can 
be understood as an interfacial pattern with a Laplacian 
instability Jll| , qualitatively similar to other Laplacian 
growth problems [T^ . For future quantitative analyti- 
cal progress, we identify the electric screening length as 
a relevant regularization mechanism. Our finding casts 
doubts on the existence of a stationary mode of streamer 
propagation with a fixed head radius. 
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FIG. 1. Evolution of spontaneous branching of anode directed streamers in a strong homogeneous background field at times 
t = 300, 365, 420 and 450. Model, initial and boundary conditions are discussed in the text. The planar cathode is located at 
2 = and the planar anode a.t z — 2000 (shown is < 2: < 1400). The radial coordinate extends from the origin up to r = 2000 
(shown is < r < 600). The thin lines denote levels of equal electron density a with increments of 0.1 or 0.2 as indicated by 
the labels. The thick lines denote the higher electron density levels 1., 2., 3., 4., 5. and 6. These high densities appear only at 
the last time step t — 450 in the core of the new branches. 
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We investigate the minimal streamer model, i.e., a 
"fluid approximation" with local field-dependent impact 
ionization reaction in a non-attaching gas like argon or 
nitrogen In detail, the dynamics is as follows: 

(i) an impact ionization reaction in local field approx- 
imation: free electrons and positive ions are generated 
by impact of accelerated electrons on neutral molecules 

drUe -I- Vr • je = drUi -|- Vr • = l^e^^el a{\£\/Eo); 

rie^i and je,i are particle densities or currents of elec- 
trons or ions, respectively, and £ is the electric field; in 
all numerical work, we use the Townsend approximation 
ao a(|5|/i?o) = OLQ exp(— £'o/|iS|) with parameters 
and Eq for the effective cross-section. 

[ii) drift and diffusion of the charged particles in the local 
electric field je — ^fieSn,; — Z?eVRne, where in anode- 
directed streamers the mobility of the ions actually can 
be neglected because it is more than two orders of mag- 
nitude smaller than the mobility /ig of the electrons, so 
j,: = 0, 

(Hi) the modification of the externally applied electric 
field through the space charges of the particles according 
to the Poisson equation Vr • £ = e{ni — ne)/eo- It is this 
coupling between space charges and electric field which 
makes the problem nonlinear. 

The natural units of the model are given by the ioniza- 
tion length Rq = ttp ^, the characteristic impact ioniza- 
tion field Eq, and the electron mobility iie determining 
the velocity vq — fJ-eEo and the time scale tq = Rq/vq. 
Hence we introduce the dimensionless coordinates 
r = R/i?o and t = t/tq, the dimensionless field E = 
£/Eq, the dimensionless electron and ion particle den- 
sities a = Ue/no and p = rii/no with no — SoEo/ieRo), 
and the dimensionless diffusion constant D — De/ {RqVq). 
After this rescaling, the model has the form: 

dt a ~ V ■ {a'E + DVa) = /(|E|) , (1) 
dtp = a/(|E|) , (2) 
p-a = V-E , E = , (3) 

/(|E|) = |E| a(|E|) ( = |E| e-i/l^l in sim.) . (4) 

In the simulations presented here, a planar cathode is 
located at z = and a planar anode aX z = 2000. The 
stationary potential difference between the electrodes 
A$ — 1000 corresponds to a uniform background field 
E = —0.5 Gz in the z direction. For nitrogen under nor- 
mal conditions with effective parameters as in , this 
corresponds to an electrode separation of « 5 mm and 
a potential difference of « 50 kV. The unit of time tq is 
K. 3 ps, and the unit of field Eo is « 200 kV/cm. We 
used D = Q.l which is appropriate for nitrogen, and we 
assumed cylindrical symmetry of the streamer. The ra- 
dial coordinate extends from the origin up to r = 2000 
to avoid lateral boundary effects on the field configura- 
tion. As initial condition, we used an electrically neutral 
Gaussian ionization seed on the cathode 



cr(r, z, t = 0) = p(r, z, t = 0) = 10"^ e 



-(z^+r^)/100^ 



(5) 



The parameters of our numerical experiment are essen- 
tially the same as in the earlier simulations of Vitello et 
al. 10 1, except that our background electric field is twice 



as high; the earlier work had 25 kV applied over a gap of 
5 mm. This corresponded to a dimensionless background 
field of 0.25, and branching was not observed. 

In Fig. 1 we show the electron density levels at four 
time steps of the evolution in the higher background field 
of 0.5. We observe that at time t = 420, the streamer 
develops instabilities at the tip. At time t = 450, these 
instabilities have grown out into separate fingers. Be- 
cause of the imposed cylindrical geometry, the further 
evolution after branching ceases to be physical. On the 
other hand, the main effect of the unphysical symmetry 
constraint is to suppress all linear instability modes that 
are not cylindrically symmetric. Hence in a fully 3D sys- 
tem, the instability will develop even earlier than here. 

Further simulations show: (a) branching does not oc- 
cur in a system of the same size in the lower background 
field of 0.25, in agreement with |l^]. (6) Branching is 
not due to the proximity of the anode, since in a sys- 
tem with twice the electrode separation (with the an- 
ode at z = 4000) and with twice the potential difference 
(A$ — 2000) — so with the same background field — , 
the streamer branches in about the same way after about 
the same time and travel distance, (c) The phenomenon 
is not specific to the particular initial condition (5). (d) 
Branching does somewhat depend on the numerical dis- 
cretization. A wider numerical mesh leads to a higher 
effective noise level; and the branching then is triggered 
somewhat earlier, (e) Occasionally, we observe a differ- 
ent tip splitting mode. In Fig. 1 at time t = 450, the 
finger on the axis develops the strongest with a exceed- 
ing 6, while in the fingers off the axis, a stays below 3. 
In the other branching mode, the first finger off the axis 
outruns the finger on the axis. 

Before we discuss the physical nature of the instabil- 
ity, we explain our numerical approach: we used uni- 
form space-time grids with a spatial mesh of 1000 x 1000. 
The spatial discretization is based on local mass bal- 
ances. The diffusive fluxes are approximated in stan- 
dard fashion with second order accuracy. For the con- 
vective fluxes a third order upwind-biased formula was 
chosen to reduce the numerical oscillations that are com- 
mon with second order central fluxes. Such oscillations 
can be completely avoided, e.g., by flux-limiting, but pre- 
liminary tests showed that the upwind-biased formula al- 
ready gives sufficient numerical monotonicity and is much 
faster. Time stepping is based on an explicit linear 2-step 
method, where at each time step the Poisson equation is 
solved by the FISHPACK routine. References for these 
procedures can be found in 
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t = 300 



t = 365 




To understand now why and at which stage the 
streamer develops a tip sphtting instabihty, in Fig. 2, 
we zoom into the streamer head. Shown are the first two 
time steps from Fig. 1 with the electron density levels 
again as thin lines, and additionally with the equipoten- 
tial lines as thick lines. One observes that during the 
temporal evolution prior to branching, both the curva- 
ture and the thickness of the ionization front decrease. 
So the width of the front becomes much smaller than 
its radius of curvature, and an interface approximation 
becomes increasingly justified. The electric field inside 
the streamer head also decreases, so that the ionization 
front more and more coincides with an equipotential sur- 
face. In summary, the ionization front evolves towards 
a weakly curved and almost equipotential moving ion- 
ization boundary. At the same time, the electric field 
immediately ahead of the streamer increases. 

We argue now that a transient stage of an ap- 
proximately equipotential and weakly curved ionization 
boundary leads to tip splitting. Conversely, we argue 
that tip splitting in the lower background field of 0.25 is 
not observed within the presently and previously [|lO| in- 
vestigated gap lengths because the transient stage of Fig. 
2 is not reached before the streamer reaches the anode. 

In fact, the streamer in Fig. 2 approaches the limit 
of "ideal conductivity" : the conducting body has <I> « 
const., while in the non-ionized region — due to 
the absence of space charges. The boundary between 
the two regions moves approximately with the drift ve- 
locity Vf (X V$ or with the diffusion corrected velocity 

Vf cx V$ (^l + 2^Da(|V$|)/|V$|) Our simu- 

lations are the first numerical evidence that the "ideal 
conductivity" limit can be approached within our model. 



This limit of ideally conducting streamers in an elec- 
tric field that becomes uniform far ahead of the front 
was studied by Lozansky and Firsov j^. They realized 
that uniformly propagating paraboloids of arbitrary ra- 
dius of curvature are solutions of this problem. They 
did not realize that these paraboloids are mathemati- 
cally equivalent to the Ivantsov paraboloids of den- 
dritic growth found earlier. The uniformly propagating 
Ivantsov paraboloids in the early 80'ies were identified as 
dynamically unstable. This is generally the case for such 
so-called Laplacian growth problems without a regular- 
ization mechanism. Since ideally conducting streamers 
also pose such a Laplacian growth problem the dy- 
namical instability of the structure shown in Fig. 2 can 
be expected, and it actually occurs as can be seen in Fig. 
1. This explains qualitatively why tip splitting occurs. 

For a quantitative analysis, a system specific regular- 
ization mechanism has to be found Its identifi- 

cation is intricate because negative streamer fronts are 
so-called pulled fronts whose dynamics is dominated by 
the leading edge rather than the nonlinear interior of the 
front Therefore standard methods like the pertu- 

bative derivation of a moving boundary approximation 
for the model does not work [|5|. (Pulling also 

implies that standard numerical methods with adaptive 
grids are inefficient.) However, the ionization front has 
two intrinsic length scales, a diffusion length and an elec- 
tric screening length. We therefore explore the approxi- 
mation oi D = 0. It is smooth for the velocity of negative 
fronts |11| and eliminates the leading edge, and hence 
suppresses the pulled nature of the front. Rather the 
front becomes a shock front for the electron density, while 
the intrinsic length scale of the electric screening layer 
behind the shock remains. 
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As a first step to understand the short wave length reg- 
ularization of perturbations due to this screening length, 
we have investigated the transversal instability modes of 
a planar ionization front in the limit _D = in a field that 
approaches the uniform limit E — —Eoo far ahead of 
the front. The planar unperturbed front propagates with 
velocity v = Eao, which equals the drift velocity of the 
electrons precisely at the shock front. The implicit ana- 
lytical front solution can be found in . In a comoving 
frame ^ — z — vt, we denote it by ( cro(Oj Po(C), ^o(0 )■ 
The Fourier components ( ak,Pk,^k ) of a transversal 
linear perturbation are defined through 



/ dkakiOe 



ikx-\-st 



+ 



etc. 



(6) 



For the derivation of the boundary conditions on the 
shock front, it is more convenient to write a single Fourier 
component as cr = cto - e''^^+'**) + crfc(^) e''^^+'** + . . .. 
With this ansatz and the auxiliary field Tpk = d(,4'ki the 
Fourier components solve the inhomogeneous equation 
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The boundary conditions at the shock ^ = can be ob- 
tained from the analytical solution in the non-ionized 
area, and from the boundedness of the charge densities: 



(8) 



The other boundary conditions are obtained by impos- 
ing that at ^ ^ — oo the electric field decays and the 
densities become constant. 

These equations together with the boundary condi- 
tions define an eigenvalue problem for s = s{k,v) with 
V — Eoo- It can be solved numerically by shooting from 
^ = towards — oo. In agreement with analytical limits 
— details will be given elsewhere — , we find 
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(9) 



This means that the electric screening length 1/Q!(|i?oo|) 
does regularize the instability of short wave length per- 
turbations from linear growth in k to the saturation value 
s{k) = \Eoo \ a{Eoc)/2. A small positive growth rate re- 
mains, but the analytical derivation of hints to the 
unconventional possibility that sufficiently curved fronts 



actually are stable to short wave length perturbations. 
This question is presently under investigation. If true, it 
would identify a most unstable wave length determining 
the width of the fingers that emerge after tip splitting. 

In conclusion, we have presented numerical evidence 
that anode-directed streamers in a sufficiently strong, but 
uniform field can branch spontaneously even in a fully de- 
terministic fluid model. We have argued that this hap- 
pens when the streamer approaches the limit of ideal con- 
ductivity. We have established a qualitative mathemat- 
ical analogy with tip splitting of viscous fingers through 
the concept of Laplacian growth, and we have analyt- 
ically demonstrated that the electric screening length 
leads to an unconventional regularization. This opens 
the way to future quantitative analytical progress. 
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